Raw data on PM2.5 has been provided by the German environmental protection agency (UBN). …
Data on total daily COVID-19 cases and deaths has been provided by the Robert-Koch-Institut. The data has been smoothed using a gaussian loess function with span 0.3 to reduce variation due to delayed reporting over the weekend and especcially on Sunday and Monday. The daily number of new infections has been computed from the daily variations of the total confirmed numbers.
The following figure shows the country wide average of daily PM2.5 and new, smoothed COVID-19 cases. Black vertical dotted lines represent the start of the contact restrictions (March 14th) and shut down (March 17th).
Some studies compare the development of atmospheric parameters and COVID-19 cases over the entire avialable time series. As both air quality and COVID-19 rates generally decrease after a shutdown event, at least parts of the identified correlations are likely caused by this external event, especially in regions with a strong influence of local activity on local air quality.
To focus on the corelated development of PM2.5 and COVID-19 cases during the early and exponential growing phase, we restrict the time series to the maximum incubation phase of 14 days prior the first reported infection and the turning point of the infection dynamics shortly after the absolute maximum. For the following figure, the time series has been restricted to this period and the daily new, smothed COVID-19 cases have been detrended using a poission regression.
Based on that, the cross-correlation between PM2.5 and detrended daily new COVID-19 cases shows that PM2.5 is both leading up to 6 days and lagging up to 9 days.
Since the dentrended time series is still rather non-stationary and to get a better idea of the time periods and date ranges related to certain time lags, a wavelet coherence analysis is performed with a loess smoother.
The analysis shows that for the period of 4 to 5 days, the PM2.5 time series is leading arround March 1 with up to 2 or 3 days (arrows towards right upward). The situation changes towards April 1st. Here the COVID-19 cases are taking over the lead, with a small time lag at the period of arround 4 days and a lag of about 2 or 3 days at a period of about 8 days.
The following figure shows clusters with similar development of daily COVID-19 cases.
The average of daily PM2.5 and smoothed new COVID-19 cases along with their detrened series within each cluster is shown in the figure below.
Based on that, the cross-correlation between PM2.5 and detrended daily new COVID-19 cases per cluster is shown below.
Finally, this is the wavelet coherence analysis per cluster.